function [channel2]=fading_corchannel(v,fc,samplingrate,num_perframe,power,gauss_seed,Tr,Rr,tx,rx,path_num,k,K,run_option,user_index)
global sim_consts;
MU_option = sim_consts.MU_option;
fdmax=v*fc/(3.6*3*10^8);
%fdmax=63.3333
%fdmax=16.2037;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
dot=2^11;
intep=samplingrate/fdmax/2-1;
p=2^8-1;
ratio=1/(intep/p);

clear intep;

number=ceil(num_perframe*ratio);

sampled_Dopler=Dopler_pdf(dot,fdmax);

[temp_noise1,temp_noise2]=gauss(dot,sampled_Dopler,k,gauss_seed);
clear sampled_Dopler;
if MU_option == 1
    if run_option == 0            
        temp1=corchanneltest3(num_perframe,temp_noise1,temp_noise2,dot,p,ratio,number,Tr,Rr,tx,rx,path_num,k,K)*sqrt(power(k)/sum(power));
        save(strcat('.\channeldata\User',num2str(user_index),'\','temp1','K',num2str(K),'k',num2str(k),'tx',num2str(tx),'rx',num2str(rx)),'temp1');
     else
        load(strcat('.\channeldata\User',num2str(user_index),'\','temp1','K',num2str(K),'k',num2str(k),'tx',num2str(tx),'rx',num2str(rx)),'temp1');         
     end;
else
    if run_option == 0            
        temp1=corchanneltest3(num_perframe,temp_noise1,temp_noise2,dot,p,ratio,number,Tr,Rr,tx,rx,path_num,k,K)*sqrt(power(k)/sum(power));
        save(strcat('./channeldata/','temp1','K',num2str(K),'k',num2str(k),'tx',num2str(tx),'rx',num2str(rx)),'temp1');
%         save(strcat('.\channeldata_epa\','temp1','K',num2str(K),'k',num2str(k),'tx',num2str(tx),'rx',num2str(rx)),'temp1');
     else
        load(strcat('.\channeldata\','temp1','K',num2str(K),'k',num2str(k),'tx',num2str(tx),'rx',num2str(rx)),'temp1'); 
%         load(strcat('.\channeldata_epa\','temp1','K',num2str(K),'k',num2str(k),'tx',num2str(tx),'rx',num2str(rx)),'temp1');  
     end;    
end;

clear temp_noise1;
clear temp_noise2;
channel2=temp1;
clear temp1;